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Abstract. The ground state of the quantum rotor model in two dimensions with 
random phase frustration is investigated. Extensive Monte Carlo simulations are 
performed on the corresponding (2+l)-dimensional classical model under the entropic 
sampling scheme. For weak quantum fluctuation, the system is found to be in a phase 
glass phase characterized by a finite compressibility and a finite value for the Edwards- 
Anderson order parameter, signifying long-ranged phase rigidity in both spatial and 
imaginary time directions. Scaling properties of the model near the transition to 
the gapped, Mott insulator state with vanishing compressibility are analyzed. At 
the quantum critical point, the dynamic exponent Zdyn — 1.17 is greater than one. 
Correlation length exponents in the spatial and imaginary time directions are given 
by 1/ ~ 0.73 and v^^ ~ 0.85, respectively, both assume values greater than 0.6723 of 
the pure case. We speculate that the phase glass phase is superconducting rather than 
metallic in the zero current limit. 
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1. Introduction 

Understanding the macroscopic state of a system of interacting bosons at zero 
temperature is of interest to the solution of a number of problems in condensed 
matter physics [U El EJ H]. The theory of superconductivity is built on Cooper pairs 
which can be treated as bosons. Mapping of flux-lines in type-II superconductors 
to a two-dimensional (2D) bosonic system has led to a deeper understanding of the 
I-V characteristics of cuprate superconductors [5], E]. More recently, the realization 
of Bose-Einstein condensation (BEC) in dilute atomic alkali gases has provided an 
experimental means to systematically explore various types of macroscopic quantum 
states and transitions between them, enriching our knowledge about equilibrium and 
dynamic properties of strongly correlated quantum systems[7l [U [9]. 

Previously, it has been established that repulsive bosons in restricted geometries 
(such as those confined by an optical lattice) may undergo a zero-temperature quantum 
phase transition from a superfluid state to a Mott insulator state as the strength of the 
interaction is increased [10]. The superfluid state is characterized by its well-known off- 
diagonal long-range order and a finite phase rigidity, whereas the Mott insulator state 
is characterized by zero compressibility and gapped particle excitations. Transition 
between the two gives over to intervening states (e.g., Bose glass or Mott glass) when 
disorder, either in the form of random on-site potential or random hopping coefficients 
are introduced [TOl [IH [IH [13] . 

In this paper we focus on a different type of disorder, i.e., random phase frustration 
on the ground state properties of a 2D interacting bosonic system. Such a situation 
has been realized experimentally in positionally disordered Josephson junction arrays, 
where the frustration can be tuned by varying the strength of a transverse magnetic 
field [31 [H]. At sufficiently strong disorder, it has been suggested that the superfluid 
state changes into a new state of matter, the "phase glass" phase [15|. Although the 
existence of such a state is not much disputed, its precise physical properties have not 
been well established [T5l [16]. 

The classical limit of the above problem can be described by a 2D XY model with 
random phase shifts[l7]. The phase diagram of the classical model has been worked out 
by Nattermann et a/. [18] using renormalization group methods. When the frustration 
is weak, only a small number of localized and tightly bound vortex-antivortex pairs 
are present in the classical ground state, so that long-ranged phase order is preserved 
at sufficiently low temperatures. However, as the frustration exceeds a certain critical 
value, the ground state becomes unstable against free vortex excitations due to a large 
distance instability [HI [iHl [20l [21]. Consequently, vortices and antivortices proliferate 
and destroy the ordered state at all temperatures. 

The gauge glass model represents an extreme case where the random frustration 
attains its maximum strength [22l [23] . Despite extensive studies by a number of 
groups [231 [211 EHl [261 [271 [281 [211 [30], the low temperature properties of the 2D model 
are still controversial. At the heart of the discussion is whether a certain form of glass 
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rigidity survives despite the presence of a finite density of free vortices in the system. 
Numerical simulations of the gauge glass model have yielded contradictory conclusions 
with regard to a finite temperature glass transition. In Ref. [31], we have carried out 
an explicit analysis of vortex configurations in the ground and low-lying excited states 
of a corresponding Coulomb gas model, where the random frustration is represented by 
a set of randomly oriented dipoles which interact with the vortices [T8| [21] . The main 
conclusion of the analysis is that, unlike the case of an Ising spin glass, the ground state of 
the 2D Coulomb gas in the random dipolar field is quite unique. The low-lying excited 
vortex states can be described in terms of a dilute gas of localized vortex-antivortex 
pairs super-imposed on a complex and critical but otherwise innocent ground state 
vortex configuration. Due to the disorder, the density of states of these pairs is finite 
at zero excitation energy. Consequently, these pairs are able to participate in dielectric 
screening under thermal equilibrium conditions. Renormalization group arguments show 
that the ground state of the gauge glass is a critical state with a phase rigidity that 
decays algebraically with distance. Thermal excitations of low energy vortex-antivortex 
pairs at temperature T lead to a finite correlation length ^ (T) which diverges as T tends 
to zero. These predictions are confirmed by direct Monte Carlo simulations of the gauge 
glass model. 

Since the power-law decay of phase rigidity with distance r requires dielectric 
screening by excited vortex-antivortex pairs of size comparable to r and the 
corresponding vortex movement at smaller scales, it is legitimate to ask if such a 
complex relaxation process can be realized in dynamic simulations and experiments. 
The slow dynamics for the creation and annhilation of distant vortex-antivortex pairs 
may indeed give rise to an apparent phase rigidity bigger than its equilibrium value. It 
is plausible that this type of glassy behavior is responsible for the previously reported 
finite temperature transition[lll [28], [29] [30], though one needs to work out the relevant 
energy and time scales for a detailed verification of this scenario. 

Quantum fluctuations of the phase lower the core energy of vortices and antivortices 
and facilitate their delocalization through quantum tunneling. Sufficiently strong 
fluctuations of this type destroy long-ranged phase ordering as in the pure case, giving 
rise to the Mott insulator state. In the following we present detailed numerical results 
to show that, when the fluctuations are weak, the ground state of the system retains 
phase rigidity both in space and in time as in the classical case. This phase glass phase 
is characterized by a flnite Edwards- Anderson order parameter, a flnite compressibility, 
but a vanishing superfluid density or helicity modulus. We also carry out a scaling 
analysis of the transition to the Mott insulator state. Our results for various critical 
exponents are signiflcantly different from those of the pure case, suggesting that the 
phase glass to the Mott insulator transition belongs to a different universality class than 
the three-dimensional (3D) XY mo del [H [2], [10!. 

The paper is organized as follows. In Sec. 2, we introduce the model Hamiltonian 
and the mapping to the (2+l)-dimensional classical model. The procedure for 
performing Monte Carlo simulations under an entropic sampling scheme is outlined. 
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Various quantities characterizing ordering in the phase glass phase are introduced. We 
also review briefly a general scaling theory developed by Fisher et al. fTO] for the analysis 
of the critical behavior at the quantum phase transition of a bosonic system. Section 3 
contains results from extensive Monte Carlo simulations of the (2+l)-dimensional model. 
The critical exponents are extracted based on a flnite-size scaling analysis. Signiflcance 
of these flndings are discussed briefly in Sec. 4. The mapping between the 2D quantum 
model and the (2+l)-dimensional classical model is presented in the Appendix for easy 
reference. 

2. The randomly frustrated quantum rotor model and its basic properties 

2.1. The Hamiltonian 

We consider a 2D Josephson junction array of superconducting grains in a transverse 
magnetic fleld. In a coarse-grained description, the Hamiltonian can be written as, 



Here Ui and 6i are the number operator of (excess) Cooper pairs and the phase of the 
superconducting order parameter on grain i, respectively. They satisfy the commutation 
relation [6j,nk] = iSj^k- In the ^-representation, we may write nj = —id/ddj. The 
flrst term on the right-hand-side of ([1]), which sums over all N sites of a square 
lattice, represents on-site Coulomb repulsion between the Cooper pairs whose strength 
is specifled by the charging energy Eq = 4e^/C, with C being the capacitance of a single 
grain. The second term represents Josephson coupling between neighboring islands with 
strength Ej, and the sum is over all nearest neighbor bonds of the lattice. 

The phase shifts in Eq. ([1]) are related to the vector potential A of an external 
magnetic fleld through. 



Here $o = he/ (2e) is the elementary flux quantum. In the present paper we shall focus 
on the maximally frustrated case where the random phase shifts uniformly 
distributed on the interval [0, 27r) and uncorrelated from bond to bond. The case aij = 
corresponds to the well-known quantum rotor model[2]. 

The Hamiltonian ([T]) is invariant under a global rotation, 6i —>■ 6i + c, all i. This 
symmetry is important in the discussion of the low energy excitations of the system. 

2.2. Mapping to a (2+l)-dimensional classical model 

It is well-known that, as far as the thermal equilibrium properties are concerned, the 
quantum rotor model can be mapped to a (2+l)-dimensional classical system using 
the Trotter formulafU [TT] . This mapping forms the basis of our quantum Monte Carlo 
calculation. 



(1) 




(2) 
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Let us start with the partition function 

Z(/5)=Tr exp{-PH), (3) 

where (3 = l/lksT) is the inverse temperature. Following the procedure described in 
Appendix A, we arrive at a classical model defined by the action flA.9|) . The model 
can be brought into a dimensionless form through the introduction oi z = t/tq as the 
coordinate in the third direction, where tq = h/ ^EqEj. After the transformation, we 
obtain 

Z{P) = J [V9] exp [--i J^"^' dzHp] , (4) 



where = /3h/To = ^EjEq/kBT, K = ^Eq/Ej, and 

= ^ E($)' - E ^os{e, - - «^.)- (5) 



2 dz ,. 

In our Monte Carlo simulations of the quantum rotor model, we choose a 
discretization along the z-axis with dz = 1 and approximate dOj/dz by 6j{z-\- 1) — 6j{z). 
This procedure appears to be rather adequate for the exploration of the phase glass 
phase and the transition to the Mott insulator state. With this choice of dz, the critical 
point Kc of the discrete model is expected to be somewhat different from that of the 
original model, though we believe the large-distance properties are not affected. The 
boundary conditions in the xy-plane are chosen to be periodic, while Eq. ( lA.lOl) is used 
along the z-direction. With these specifications, we obtain a classical (2+l)-dimensional 
lattice model where K plays the role of a "quantum temperature" . As such, the model 
can be simulated using the entropic sampling scheme [32]. Since the algorithm allows the 
system to explore configurations over a broad range of values of the classical action ([5]) in 
a single simulation, it has a better chance to generate statistically independent samples 
for the calculation of thermal averages, which may become a problem for conventional 
Monte Carlo methods in the glass phase. In addition, once equilibrated, values of any 
measurable quantity over a broad range of K values can be readily calculated. This is 
particularly useful for analyzing the zero-temperature quantum phase transition. More 
details about our implementation of the scheme can be found in Ref. [33j . 

2.3. Helicity modulus and compressibility 

Helicity modulus can be related to the superfiuid density and hence provides a direct 
measure of the superconducting order. It is defined by considering the change of the 
free energy under a phase twist across the system. Quite generally, such a twist modifies 
the action to Hp = Hp + 6Hp. The change in the free energy, to the second order in 
6Hp, is given by, 

5F= -kBT\n{Z/Z) 

= — fc^T ln(exp(— — / dz6Hp)\ 
\ K Jo I 
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Here (•) denotes thermal average. For a twist A^. per bond in the x-direction, we have 

6Hp = A^Yl sin(6'i+^ - Oi + Oj,^.) + ^A^ ^ cos(6'i - 6*,;+^ - ai^r,) + ... (7) 
i ^ i 

Writing SF = ^NpsEjAl, we obtain the following expression for the helicity modulus, 

N r^'^ - - - o 

Ps = {cosie, - ei+, - a,,,.)) - ^ / dz[{.uo)Uz)) - (JM) ]- (s) 

K Jo 

Here Jx{z) = N^'^ J2i s^^{(^i+x{z)—9i{z)+ai^x) is the average current along the x-direction 
in layer z, and the overline bar denotes spatial average. 

The compressibility k is related to the change in free energy for a phase twist A in 
the z-direction|TO]. The corresponding change in Hp is given by, 



SH^^AE'^.InA^ (9) 

Inserting Eq. into Eq. ([6]), and noting that Jq"" dz{d6i/dz) = 6i{Lz) — 9i{0) = 27mi, 
where rij is the number of turns the angle on site i makes along the z-direction, we 
obtain 

SF . InE,^' - ^{(l.p.n,f) = iiV£,KA^ (10) 

where 

Here (ui) =0 due to symmetry. 

2.4- The Edwards-Anderson order parameter 

The Edwards-Anderson order parameter for the quantum rotor model can be defined 
via 



qea = lim (e^[^.(*)-e.(o)]), (12) 

t — ^oo 

where e*^-'*-*^ = ^-^Ht/n^iej ^iut/h 

Consider now the auto-correlation function in imaginery time. 



^^y^-{l3h-T)H/h^iej^-HT/h^- 

Z 



where z = t/tq, and the last average is carried out in the (2+l)-dimensional classical 
ensemble. At T = 0, the long-time limit t ^ oo can be replaced by the limit r — oo. 
Hence we may write, 

gEA= lim C{LJ2,L,). (14) 

L,— >oo 
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2. 5. Scaling properties near the Mott insulator transition 

The Mott insulator state corresponds to the "high temperature" phase of the (2+1)- 
dimensional classical model where quantum phase fluctuations lead to a vanishing 
superfiuid density, vanishing compressibility, and an exponentially decaying phase- 
correlation function in imaginary time. Due to the anisotropic form of ([5]), two different 
lengths ^ and are needed to describe the spatial and temporal correlations in the 
system. As K approaches its critical value Kc at the transition from the Mott insulator 
to the phase glass, both quantities are expected to diverge as 

^^IK-K.r, ^ \K - K,r% (15) 

where and z/^ are the respective exponents. Since corresponds to a characteristic 
energy or frequency scale in the quantum rotor model, the ratio z^yn = ^z/^ defines 
the dynamical exponent at the transition. In the phase glass phase, both ^ and are 
expected to be infinite. The exponent z^^^ of the phase glass phase, which may be 
different from its value at the transition, can be determined from suitable finite-size 
scaling properties. 

Earlier, Fisher et a/. [10] proposed a general scaling theory for the quantum phase 
transition in bosonic systems. In particular, based on the assumption that the singular 
part of the free energy (or ground state energy in the quantum model) in a correlated 
volume ^'^^z is of order i^c? they determined the scaling dimensions of the superfiuid 
density and compressibility in the transition region, 

K^\K-K,\''% (16) 

where d is the spatial dimension, and 

C = {d-2)v + Vz, Cz = du- Vz (17) 

are the scaling exponents. 

The auto-correlation function (|T3l) . on the other hand, is expected to decay as a 
power law at K = K^^ 

C(r, oo) ~ ^-i-(rf-2+r?)/^dyn^ (18) 
where is a new exponent. Scaling arguments then yield, 

gEA ~ |^-^c|^'^~'+''^""^"^ (19) 

For the (2+l)-dimensional classical model, the "specific heat" cyiK^ = 
-Kd'^f/dK^ with / = F/{NLz) is expected to exhibit singular behavior at Kc. From 
the above assumption for the singular part of the free energy ~ one obtains 

the specific heat exponent, 

a = 2 - du - Uz- (20) 

In the pure case = and = 2, the quantum rotor model is mapped to the 
3D classical XY model. Consequently 2;dyn = 1, and C = Cz = = '^z, i-e., ps and 
K scale the same way as E,^^ which provides the only energy scale of the problem. 
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Numerical calculations have yielded t'soxY — 0.6723[31]. From the scaling relation 
7 = (2 — ri)^ ^ 1.319 we obtain t^sdxy — 0.04. The exponent asDXY — —0.017 is 
also very small. One of our numerical tasks below is to check whether the same set of 
exponent values apply to the phase glass to the Mott insulator transition. 

3. Simulation results 

We have carried out extensive Monte Carlo simulations of the (2+l)-dimensional 
classical model ([5]) under the entropic sampling scheme. The system is chosen to be 
a cubic lattice of Lz layers each containing N = sites in the x|/-plane. Test runs were 
performed on the unfrustrated quantum rotor model (referrd to below as the pure model) 
which generated results in good agreement with previous studies on the superfluid to 
the Mott insulator transition at i^co — 2.55. Unless explicitly stated, the data presented 
below for the disordered case are obtained from averages over 30 to 100 samples at any 
given size. This seems to be sufficient for illustrating the behavior of the phase glass 
phase and for determining the critical exponents of the transition within the limit set 
by the system size we are able to investigate. 

Figure 1(a) shows the "specific heat" data for the pure model against K for six 
different system sizes. The cusp singularity at K^q with a very small exponent a is 
evident. In comparison, as seen in Fig. 1(b), the singularity is much weaker in the 
randomly frustrated model. The dashed line in the figure with 

a ~ -0.3 (21) 

indicates a possible behavior in the infinite size limit that is consistent with our data, 
though in general the critical amplitudes on the two sides of the transition need not be 
the same. 

Figure 2 shows the helicity modulus ps and the compressibility k of a 16'^ system 
against K for four different disorder realizations. Both quantities become vanishingly 
small when K > Kc — 1.98. At smaller values of K, ps is strongly sample-dependent 
and takes on both positive and negative values. For a given sample, ps may also be a 
non-monotonic function of K. The disorder-averaged value of ps, on the other hand, 
remains close to zero. A nonzero ps for individual samples signifies "freezing" of vortex 
loops that enables long-ranged phase rigidity to develop. To understand the origin of 
negative values for ps, one may consider a more general twist boundary condition in 
each layer, O^+i^y = O^^y + and O^^y+L = Gx,y + Ay (see, e.g., Ref . [25] ) . Due to the 
random phase shifts, the minimum of the free energy is in general achieved at some 
nonzero values of and Aj,. Depending on the particular choice of the random phase 
shifts, the curvature ps of the free energy at A^^ = A^ = may take on positive or 
negative values. A sign change can occur when one or more vortices are relocated on 
the scale of the system size. Such a move does not alter the local phase gradients in a 
significant way and hence costs only a small amount of energy. As K is varied, one may 
envisage a change in the relative statistical significance of configurations that differ in 
this way, leading to the observed non-monotonic behavior. 
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Figure 1. The "heat capacity" against K = y/ Eq /Ej for (a) The pure model and (b) 
the randomly frustrated model at six different system sizes as indicated in the figure. 
The dashed line in (b) represents a possible behavior at infinite system size. Here 
Kr = 1.98. 




In contrast, the compressibilities of the four samples differ only slightly from 
one another, with no significant broadening near K^. Other quantities, such as the 
specific heat and correlation functions, show the same behavior. Therefore the system 
is expected to be self-averaging in the large-size limit. 



Phase glass and zero-temperature phase transition. 



10 



The critical exponents for the transition can be determined by applying suitable 
finite-size scaling forms. Below we shall perform the analysis using as the scaling 
variable. The general finite-size scaling ansatz of a quantity X then reads, 

X{K, L, L,) = L-^I^^X{{K - K,)Ll/''%L''^y-/L,), (22) 

where x is the critical exponent for the quantity X. Equation (!22|) is quite difficult to 
use in general due to the simultaneous presence of two scaled variables. However, as we 
shall see below, the exponent Zdyn is quite close to (but larger than) one so that, for the 
range of system sizes considered, the second argument is approximately constant and 
does not affect significantly the value of the function. The validity of this assumption 
is justified a posteriori by the consistency of the exponent values obtained. 

Figure 3(a) shows the disorder averaged compressibility k against K for six different 
system sizes, ranging from L = = 4 to L = = 16. At K = (as indicated by 
the dashed line in the figure), the decay of n against Lz can be fitted to a power law 
with an exponent (z/'^z = 0.7 ± 0.1. From the scaling relation (fT7|) we obtain 

2dyn = 1.17 ±0.07. (23) 

Combining this result with our previous estimate a = 2 — 2u — Uz — —0.3 yields 

u ~ 0.73, Uz ~ 0.85, Cz = 2u-Uz- 0.61. (24) 

The scaling plot shown in Fig. 3(b) is generated using these exponent values. A 
reasonable data collapse is seen. Note that there is a slight increase in the slope of the 
curves with increasing L or Lz, which is consistent with the expected trend associated 
with a gradual increase of the scaled variable L^'^y^/Lz for this set of data. 

We now examine the behavior of the Edwards-Anderson order parameter = 
C{Lz/2, Lz). As shown in Fig. 4(a), the effect of finite size is more pronounced here as 
compared to k in Fig. 3(a). Also, the dependence of ^ea on Lz at a finite L is complicated 
by the fact that, in a finite system, the spectrum of H is always discrete so that, strictly 
speaking, ^ea as defined by iHM vanishes in the limit (3 ^ oo. Nevertheless, for the 
choice Lz = L, we observe that data at a given K < can be well fitted by a quadratic 
function in l/L^, as shown in Fig. 4(b). Using such an extrapolation procedure, we 
obtain a value for ^ea at each K in the infinite size limit as indicated by the dashed 
line in Fig. 4(a). The approach of ^ea to zero as K tends to Kc can be described by a 
power-law with an exponent 0.8. Comparing with Eq. (fTOjl and the values mentioned 
above for u and Vz-, this result is consistent with a small value for the exponent rj. 

4. Summary and discussions 

The main findings of our numerical investigation of the quantum rotor model at maximal 
random frustration are summarized as follows. When the charging energy Eq (or boson 
repulsion) is small compared to the Josephson energy Ej (or boson hopping energy), the 
system behaves quite similarly to the classical gauge glass model at zero temperature. 
The compressibility k is positive and increases with decreasing K = JEq/Ej. The 




Figure 3. (a) The disorder averaged compressibility for six different system sizes 
L = = 4, 6, 8, 10, 12, and 16. (b) A scaling plot using the exponents discussed in 
the text. 




2.4 2.8 




Figure 4. (a) The disorder averaged Edwards-Anderson order parameter for six 
different system sizes L = Lz = 4,6,8,10,12, and 16. Dashed line indicates 
extrapolated value using a quadratic fit. (b) The dependence of qea on for selected 
values of K . 
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Edwards-Anderson order parameter gsA also assumes a finite value, signifying long- 
ranged phase ordering in time. The 1/r (or l/L^) correction to ^ea in imaginary time 
[see Fig. 4(b)] is the same as in the isotropic 3D XY model, suggesting a dynamic 
exponent = 1 and linear dispersion for the gapless phase modes. These properties 
do not seem to be affected by the power-law decaying spin-glass stiffness in the spatial 
direction obtained previously [25| [261 13T] . Our results suggest that quantum tunneling 
between the classically distinct low energy vortex states as discussed in Ref. |3T] is 
suppressed at large distances in the phase glass. The existence of spatial (glassy) phase 
order is also supported by a nonvanishing helicity modulus ps below Kc in individual 
samples, though due to the random phase shifts, ps does not have a definitive sign. 
These observations suggest that diffusive transport of vortices under an applied current 
is unlikely in the present model except perhaps at K = Kc, casting doubt on the link 
between the phase glass and bose metal when the only quenched disorder in the system 
is in the form of random phase frustrations ^15j. 

We have also attempted to determine the critical exponents characterizing the 
phase glass to the Mott insulator transition. Due to the relatively small system sizes 
available, the estimated values of the exponents should be considered as only tentative. 
With this caveat in mind, our data are broadly consistent with a dynamic exponent 
Zdyn — 1-17 ± 0.07 greater than one, and correlation length exponents u ^ 0.73 and 
Uz — 0.85, both greater than their value 0.6723 in the pure case. It would be interesting 
to confirm (or disprove) these results with more efficient numerical algorithms applied 
to the model. 
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Appendix A. Path integral representation of the partition function 



Following the standard procedure, we define M = exp{—PH/n) and write 

Z{f3) = Tr(M"). (A.l) 
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When the integer n is large, we may write 

M = exp{-l3{HQ + Hj)/n) ~ exp(-/5^Q/n) exp(-/?^j)/n), 

where 
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(A.2) 

(A.3) 
(A.4) 



Hj = - EjY^ cos{9i - 9j - aij), 

are the Coulomb and Josephson energies, respectively, which do not commute. 

Let I {6'}) be a state where each rotor j has a definitive phase 6j G [0, 27r), and 
be a state where each rotor j has a definitive angular momentum rrij = 0, ±1, ±2, . . . 
The matrix elements of M can be written as, 

M{{e'},{e})^{{e'}\M\{e}) 

^ ({^'}| expi-f3HQ/n) exp{-f3Hj)/n)\{e}) 

= {{9'}\expi-PHQ/n) ^ \{m}){{m}\exp{-PHj)/n)\{9}) 



exp 



PEj 
n 



PEq 



{m} 

k - 9j - ttij) 
Y,m] + iY,mj5ej 



X E [-^ E ^? + ^ E , (A.5) 

{m} j j 

where 69 j = 9'- — 9j and we have used = exp(z mj^'j). With the help of 
the Poisson summation formula, 



E /(^) = E 

m=0,±l,±2,... 



dufiu) exp{2mus) , 



(A.6) 



s=0,±l,±2,, 

the sum over the angular momentum eigenstates {m} can be carried out, 

(3Eq 

/3Eq 

3 



E^xp 

{m} 



2^ Y.^]^^Y.^3^^3 



\ \du\ exp 

W 

A ^ exp 



2n 



n 



^(5% + 27rs 



(A.7) 



2PEq^ 

where A is a numerical constant. 

With these preparations we may carry out the matrix multiplication and trace in 
Eq. f[0|) to obtain. 



/n— 1 
[V9]l[Mmrk^i)},{0in)}), 

k=0 

1 rf^f>- 



(A.8) 
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where = khf3/n is the imaginery time coordinate and 



Hc = 7^T.^'-EjY: cos{9, - 9, - a,,) (A.9) 
^^Q j (ij) 



is the classical action. The integration over the Ojir^ys are on the infinite domain 
(— CX3, oo) for all k except at = 0, where the domain [0, 27r) is taken instead. This 
procedure takes care of the sum over the s,'s as in Eq. f[0) . Note that 

Ojiph) = ej{0) mod 27r. (A. 10) 
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